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For sparse materials like graphitic systems and carbon nanotubes the standard density functional 
theory (DFT) faces significant problems because it cannot accurately describe the van der Waals 
interactions that are essential to the carbon-nanostructure materials behavior. While standard im- 
plementations of DFT can describe the strong chemical binding within an isolated, single-walled 
carbon nanotube, a new and extended DFT implementation is needed to describe the binding be- 
tween nanotubes. We here provide the first steps to such an extension for parallel and concentric 
nanotubes through an electron-density based description of the materials coupling to the electrody- 
namical field. We thus find a consistent description of the (fully screened) van der Waals interactions 
that bind the nanotubes across the low-electron-density voids between the nanotubes, in bundles 
and as multiwalled tubes. 



I. INTRODUCTION 

Soft and sparse materials constitute an important challenge for our first-principle quantum- 
physical account of structure and dynamics. These structures are defined by voids of very-low- 
electron density regions across which material binding is communicated exclusively via the coupling 
to the electrodynamic fields, the so-called dispersion or van der Waals (vdW) forces. Traditional 
implementations of density functional theory (DFT) cannot account for this type of interaction 
because the density functional (DP) are defined as local or semilocal and fail to describe the 
truly nonlocal nature of the electrodynamic coupling. Nevertheless, the traditional approaches 
are imperative for accurate predictions of the binding within the dense regions of finite electron 
densities. The many opportunities in biological systems, in organic-molecular liquids and in the 
carbon nanostructures motivate a continuous search for a consistent combination of traditional 
DPT and vdW corrections. This biophysics and nanotechnology perspective motivated our recent 
proposal for a vdW-DP for layered systems J^, 2., J., 4] . 

The interaction in nanotube bundles, ropes |5|, and within multiwalled nanotubes as 
addressed here, represents another type of sparse-matter problem with a nanotechnological rele- 
vance 0- Nanotube bundles are often directly produced and/or extracted to align the nanotube for 
Raman spectroscopy whereas multiwalled nanotubes are produced in many fabrication techniques. 
Both the nanotube bundles and the concentric multiwalled nanotubes are stabilized by a competition 
between the kinetic-energy repulsion and the long-ranged van der Waals binding. At the relevant 
binding distances (3 to 3.5 A or 6 to 7 Bohr radii) the traditional first-principle DPT calculations 
cannot account for the combined interaction. The issues and challenges are similar to and closely 
related to the corresponding problem of ab-initio calculations of the binding in graphite Q . 

An accurate quantum-physics account of the binding in nanotube bundles and in multiwalled 
nanotubes must include a calculation of the dispersion forces defined by the electrodynamics cou- 
pling. To illustrate the importance of these van der Waals contributions we here use a simple model 
description of the single-tube electrodynamics response fitted to DFT calculation of the nanotube 
electron system and response. While simple, our approach improves the traditional Hamaker-type 
of estimates in several ways: (i) it offers a method for a first-principle calculation of the inter- 
action strength, (ii) it includes the important local-field screening effects as directly defined by the 
anisotropic spatial variation in the nanotube electron density (which we here obtain from DPT); 
and (iii) it consistently accounts for the net anisotropic dielectric response that describes extended 
molecules. In this manner we establish a link between the microscopic description (based on accurate 
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FIG. 1: Sketch of parallel nanotubes. Top panel: Tubes with radii Ri and R2, seen along the z\ and 22 
axes. Bottom panel: A closer look at the closest interaction region. Whereas all parts of a nanotube in 
principle interact with all parts of the other nanotube the strongly attractive character of the dipole-dipole 
interaction causes most of the contribution to the collective interaction to come from within the region near 
the line of closest contact. 



DFT calculations of the electron density variations and static susceptibility 3 ) and the phenomeno- 
logic descriptions, and we provide a method to test the limit of applicability. In a wider sense our 
nanotube study also provides an indirect insight on how the electron-density variation determines 
the optical response of large extended (organic) molecules as used, for example, in liquid-crystal 
displays 0. 

In the present analysis of the nanotube van der Waals interaction we focus on interaction 
str eng th at intermediate separations when the interaction is defined by the dipolar contribution 
p. llCl||. This emphasis permits a number of very important simplifications. With a set of approx- 
imations we succeed in using the high degree of symmetry which exists in the case of parallel and 
concentric nanotubes to split the interaction into an effective coupling constant (given by one fre- 
quency integration) and a nanotube geometry factor. In the worst case the computation is thus 
reduced to just a two-dimensional spatial integration, as shown in the following. Moreover, both 
for the cases of concentric nanotubes and of identical parallel nanotubes we further succeed in ex- 
pressing this general integral fully or partially in terms of spatial function for additional efhciency 
gain. 



II. THE NANOTUBE ELECTRODYNAMIC RESPONSE 



Our calculations of the physical intcrtube attraction are based on accurate first-principle DFT 
calculations of the nanotube electron density, the overall, nanotube electrostatic response, and an 
approximate treatment of the nanotube electrodynamical response, local-field effects, and resulting 
intertube attraction. The approach is described in Ref. and will be summarized here. 

Adapting the plasmon-pole model of Ref. [l^ , we approximate the local microscopic electron 
response (for a sub-atomic-scale element of the nanotube wall) by the bare dynamic susceptibility 
Xo (^(r), M, uq) = n{r)/{u'^ + ^^o) (complex) frequency u. Hartree atomic units are used with 
symbols oq for the Bohr radius and Ha for the hartree (« 27.21 eV). The nanotube electron density 
n(r) (at position r) is determined directly from first-principle DFT calculations, and the value of 
the effective frequency cut-off uq is obtained by comparing first-principle calculations of the static 
nanotube susceptibility to a calculation including the local-field effects produced by the model of 
Xo- 

From this model of xo obtain a corresponding effective susceptibility tensor XoS which 
describes the ratio of the locally induced polarization to the externally applied electric field. We 
describe the response of a nanotube to an applied electric field using a local cylindrical coordinate 
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FIG. 2: Van der Waals interaction of pairs of parallel nanotubes with pairwise identical radii, as a function 
of smallest separation A — d~2R. The attraction of large nanotubes (radius 3i?(4_4)) is greater than that of 
small tubes (radius ii(4,4)) at the same A. Large nanotubes have less curvature and thus at a given distance 
A a pair of large tubes has a larger amount of charge at distance approximately A, yielding a larger total 
binding per tube length. 




system (s, 6, z) with the z-axis along the nanotube axis and s a vector perpendicular to the z-axis. 
Approximating the nanotube electron density n(r) by its radially averaged value n(s) the local 
effective susceptibility is given by the relation 

Xcff [n{s)] (u)Eappiiod = -Xo {n{s),u, Uq) V(/)(s, u) (1) 

where the local electric field — V0(s,m) is given by charge conservation V • {(1 + 47rxo) V(/)} = 0. 
In general, the effective susceptibility tensor XoR will be anisotropic with off-diagonal elements. 
However, when transforming the tensor XcS into cylindrical coordinates all off-diagonal elements 
vanish and only the diagonal elements xfff ; P = s,9, z are left. 

For an external electric field of strength Eq(u), oriented perpendicular to the nanotubes or 
along the nanotubes, we may factorize the electric potential as 0(s, u, uq) — —Eo{u)W{s, u, uq) cos 9 
and solve V • {(1 -I- 47rxo) V^} = 0. The components of the effective susceptibility then become 

Xoff = XodsW, Xcff = XoW/s, and xk = ^o- . . 

To proceed with an efficient determination of the internanotube van der Waals interactions 
we exploit the special form of the nanotube electron density. In the thin-wall approximation we 
approximate the radial behavior of the effective susceptibility x2f[{s, u, uq) x'iR^ u, uo)S{s — R), 
by a weighted radial delta function such that / (Ps xJsi^i mq) = 2ttRx^ {R, u, mq). 

These approximations are the basis for derivation of the nanotube-nanotube interaction. 



III. THE NANOTUBE-NANOTUBE VAN DER WAALS INTERACTION 

For a pair of parallel nanotubes at center-to-center separation d (Fig. ^ we define two sets of 
local cylindrical coordinate systems with origos separated by d, and with indices 1 and 2 referring 
to the two nanotubes and their local coordinate systems. Within the dipole-dipole approximation 
we then find the van der Waals energy 

1 r°° 

^vdw - -TT / duTTace{xiTi2id)x2T2iid)} (2) 
Jo 

= E il dux^{R^,u,u^^^)x'{R2,u,u'i'^)R,R^ 

X d9i d92l^^{9i,92,Ri,R2) (3) 
Jo Jo 
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FIG. 3: Van der Waals interaction of pairs of parallel nanotubes as a function of center-of-tube to center- 
of-tube distance d. Gray lines are the asymptotes E^^-y^/L ~ —Cs/dP with Cs — Qir^ Jtot{R, R,uo,uo)/64: 
for R = -R(4,4) and R = 3-R(4^4) . 



with the geometry terms 



1 pi^/^ poo 
^ J -L/2 J -00 



(4) 



and the dipole coupling T12 — — V1V2 |r2 — riP"'^ between elements of the two nanotubes. We made 
use of the thin-wall approximations to carry out the si and S2 integrations. The geometry terms 
Ip^ are of a form for which the Z21 integration can be easily be carried out [l^ (3.241.4)]: 



dz' 



21 



"-21 



{^21 + -^21)" 



24 



'21 



(lO-^i) 



(5) 



with = 1, 3, or 5, and S21 = |s2 — Si| 



The cut-off frequencies and u\^' for the two nanotubes depend in principle on the nanotube 



radius, but as found in Ref. p4| the macroscopic susceptibility is not very sensitive to the value uq. 



,(2 



In practice we thus take the large-nanotube value of uq for all nanotubes, Uq 



(1) 



,(2) 



0.30 Ha. 



For parallel nanotubes we can thus reduce the original simultaneous six-dimensional spatial 
(ri, r2) and one-dimensional {u) frequency integral to a sum of terms each factorized into a one- 
dimensional frequency integral times a two-dimensional 62) spatial integral. The frequency 
integrals, now decoupled from the spatial integrals, do not depend on the nanotube separation d 
and are therefore identical to the frequency factors defined in Ref. p4| 



Jp^{Ri,R2,ui^\u^Q^) 



dux^iRi,u,u^^^^)x-'{R2,u,u'i^)RiR2 



(6) 



As shown earlier Tlj, the special case of concentric nanotubes (d — 0) further reduces the number 
of spatial dimensions to be integrated. 

In general it is not possible to carry out all of the nine spatial integrals If3-y{0i,62, si, S2) to 
express them in terms of known functions. The complications, however, are reduced for certain 
special cases such as concentric nanotubes and, as shown below, for parallel nanotubes of equal 
radius. 

To emphasize the effect of the cylindrical symmetry on the van der Waals interaction (O we 
would like to pursue analytic calculations as far as possible. For this reason we shall assume that the 
effective radial and tangential susceptibility are identical, = . The error of this approximation 
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FIG. 4: The geometry part of the van der Waals interaction of concentric nanotubes for fixed inner-tube 
radius -Rinncr and as a function of outer-tube radius R. 



is largest for small frequencies u, and is at most 10%. Thus, with Jgg = Jbs = JsB = Jee, Jsz = Jbz 
and Jzs — JzB, the van der Waals energy per nanotube length becomes 



vdW 
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In Eq. ITJ each sum of /-terms has the same functional form 

9tt 



d02{Jss {Iss + Is0 + les + he) 

- Jzs {Jzs + Ize) + Jzzizz} ■ 



2, Ri, R2) = (iR2Cos02 + d - Ricoseif + {R2sm02 - Risin0if 



-5/2 



(7) 

(8) 



times an integer factor: Iss + he + l0s+ lee = 3/tot, Isz +Iez = 5/tot, Izs + Ize = 5/tot, Izz = 19/tot. 
This is an important observation, as it enables the factorization into frequency and spatial integrals 
not only for each term in the sum of l(7|), but indeed for the full sum of (|7||. If we therefore define an 
effective frequency integral Jtot = iJzz + ^Jzs + 5 Jsz + ISJ^^ the van der Waals energy per length 
can be written as the product 



r — -7^Jtoti-Kl,^2,Wo '"0 ) 



d0i 



d02li 



, 6*2, i?2) 



(9) 



Eq. is our result for arbitrary-size parallel nanotubes, including concentric nanotubes. The 
two- and one-dimensional spatial and frequency integrals involved in are already a considerable 
simplification compared to the original 6-1-1 dimensional integral, and can be numerically solved. 

For identical nanotubes, Ri = R2 = R, the two spatial integrals in @ further reduce to a 
single integral over the variable £, — {02 — Si)/2 
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where 2F1 is a hypergeometric function. Replacing one integral in this way with a well-known, 
tabulated special function speeds up the numerical evaluation of E^dw and enhances the accuracy. 

In Figure [3 we plot i?vdw as a function of the smallest distance between elements of the 
nanotubes, A = d — Ri — R2 = d — 2R for three different radii. In this and the following figures we 
have assumed a constant value of uq = 0.30 Ha. 

In the asymptotic limit, d » Ri, 
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we recover to lowest order the traditional London-theory macroscopic dr^ behavior J^J of two thin, 
parallel cylinders, albeit with different higher-order correction terms because nanotubes are hollow, 
not solid, cylinders. The limit of the expansion is illustrated by the plots in Figure 13 for the 
identical tubes presented in Figure |21 

Fi gur e El shows our results for another important special case, namely the concentric nanotube 
system When the two local nanotube coordinate systems coincide {d = 0) both of the two 

spatial integrals in jnj can be solved exactly [T^ . (2.584.58)] to give -Evdw in terms of the first and 
second complete Legendre elliptic integrals K{k) and E{k) 

EvdW _ StT Jtot(-Rl, ^2, Wq^\ Uq^^) 



32 (i?2+i?l)3(i?2 



4-|i±|LiJ ( ^y^^] K ( l^^^l } . (12) 



{R2-R1? \Ri + R2j \Ri + R2 

In Figure^the geometry part of ifT^ . £'vdw/(>^toti), is shown. 

We stress that the repulsive nanotube-nanotube interaction, resulting from the kinetic-energy 
repulsion, is not treated within our formalism. 



IV. CONCLUSIONS 



In summary, we report a quantum-physics calculation of the vdW binding in nanotube bundles 
and in concentric (multiwalled) nanotubes which is based directly on first-principle calculations of 
the electron density and on the electron density response. While present DFT implementations 
cannot account for the intertube interactions we have thus identified a method for a consistent 
combination with ab-initio calculations. Together with our recent progress for the layered systems 
our results indicate that a combined vdW-DF is feasible and hence promises an integration for a 
full quantum-physics account of the sparse and soft materials used for example in biophysics and 
nanotechnology. 
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